A model for cross-cultural reciprocal interactions through mass media 
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We investigate the problem of cross-cultural interactions through mass media in a model where 
two populations of social agents, each with its own internal dynamics, get information about each 
other through reciprocal global interactions. As the agent dynamics, we employ Axelrod's model 
for social influence. The global interaction fields correspond to the statistical mode of the states of 
the agents and represent mass media messages on the cultural trend originating in each population. 
Several phases are found in the collective behavior of either population depending on parameter 
values: two homogeneous phases, one having the state of the global field acting on that population, 
and the other consisting of a state different from that reached by the applied global field; and a 
disordered phase. In addition, the system displays nontrivial effects: (i) the emergence of a largest 
minority group of appreciable size sharing a state different from that of the applied global field; (ii) 
the appearance of localized ordered states for some values of parameters when the entire system is 
observed, consisting of one population in a homogeneous state and the other in a disordered state. 
This last situation can be considered as a social analogue to a chimera state arising in globally 
coupled populations of oscillators. 

PACS numbers: 89.75.Fb, 87.23.Gc, 05.50.+q 



I. INTRODUCTION. 

The study of cross-cultural experiences through mass- 
mediated contact is a topic of much interest in the Social 
Sciences [Tj-ll] . Many of those studies have focused on the 
effects of cultural product consumption on audience be- 
liefs, emotions, and attitudes toward the group originat- 
ing these cultural products. For instance, several works 
have investigated the process by which international au- 
diences develop American values, norms and stereotypes 
about America through the experience of watching Amer- 
ican television series 0-0] • Other works have explored 
the political impact of international television across bor- 
ders Q . The expansion of broadcasting and telecommu- 
nication industries in recent times has led to an increase 
in the exchange of mass media products across countries 
and social groups. As a consequence, peoples of differ- 
ent groups that may have had little direct contact with 
each other can, however, have access to their reciprocal 
mass media messages. For example, the growth of media 
channels in East Asia has brought changing patterns of 
cultural consumption: younger generations in China are 
drawn to Korean pop stars; Korean people have begun 
to collect Chinese films; Japanese audiences await the 
broadcast of non-Japanese Asian dramas [|[. 

In the current research in complex systems, there is 
also much interest in the investigation of models of social 
dynamics Q. Many of these systems have provided sce- 
narios for investigating new forms of interactions and for 
studying new collective phenomena in non-equilibrium 



In this context, the model introduced 
17] to investigate the dissemination of cul- 



systems 
by Axelrod 

ture among interacting agents in a society has attracted 
much attention from physicists fl8l - l29l ]. In this model, 



the agent-agent interaction rule is such that no inter- 
action exists for some relative values characterizing the 
states of the agents that compose the system. This type 
of interaction is common in social and biological systems 
where there is often some bound or restriction for the 
occurrence of interaction between agents, such as a sim- 
ilarity condition for the state variable [30| - |34| . 

In this paper we investigate the problem of cross- 
cultural interactions through mass media in a model 
where two separated social groups, each with its own in- 
ternal dynamics, get information about each other solely 
through reciprocal global interactions. We address the 
question of whether two societies subject to reciprocal 
mass media interactions become more similar to each 
other or if they can mantain some diversity. Specifically, 
our system consists of two populations of social agents 
whose dynamics is described by Axelrod's model, mu- 
tually coupled through global interactions. The global 
interactions act as fields that can be interpreted as mass 
media or cultural trends [24] . [35j originated in each pop- 
ulation. 

The existence of non-interacting states in the dynam- 
ics, as well as the competition between the time scales 
of local agent-agent interactions and the responses of the 
endogenous global fields, lead to nontrivial collective be- 
haviors, such as the emergence of a largest minority group 
in a population, sharing a state different from that of the 
applied global field, and the occurrence of localized or- 
dered states. In this last case, one population reaches 
a homogeneous state while several states coexist on the 
other. This situation can be considered as a social ana- 
logue to a chimera state arising in globally coupled pop- 
ulations of oscillators .•'{<> Hj. 



In Sec. II we present the model for two interacting 
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populations of social agents and characterize the collec- 
tive behavior on the space of parameters of the system. 
The nature of the observed localized ordered states is in- 
vestigated in Sec. III. Section IV contains the conclusions 
of this work. 



II. THE MODEL. 

We consider a system of N agents consisting of two 
populations or subsets: a and /3, with sizes N a and Np, 
such that N = N a + Np. The fraction of agents in subset 
a is N a /N and that in subset (3 is Np/N. 

Each subset consists of a fully connected network, i. e., 
every agent can interact with any other within a subset. 
We employ the notation [z] to indicate "or z" . The state 
of agent i G a[/3] is given by an F-component vector 
x^tgAi), (/ = 1,2, ... ,F), where each component can 
take any of q different values x^Ji) G {0, 1, ...,q — 1}. 

Let us denote by M a = (M*, M{, M%) and 
Mp = (Mp,...,Mp,..., Mp) the global fields defined as 
the statistical modes of the states in the subsets a and 
P, respectively, at a given time. This means that the 
component -^„rgi is assigned the most abundant value 
exhibited by the /th component of all the state vectors 
x^ a rpAi) in the subset a[(3}. If the maximally abundant 
value is not unique, one of the possibilities is chosen at 
random with equal probability. In the context of social 
dynamics, these global fields can be interpreted as mass 
media messages about "trends" originated in each popu- 
lation. 

Each agent in subset a is subject to the influence of 
the global field Mp, and each agent in subset /3 is subject 
to the influence of the global field M a . Figure Q] shows 
the configuration of the two populations subject to the 
influence of their reciprocal global fields. 




FIG. 1: Representation of two populations a and j3 interacting 
through their reciprocal global fields M a and Mp, each acting 
with intensity B. 

Starting from random initial conditions in each subset, 



at any given time, a randomly selected agent in subset 
a[/3] can interact either with the global field Mpi a i or 
with any other agent belonging to a[0\. The interaction 
in each case takes place according to the dynamics of 
Axelrod's cultural model. 

The dynamics of the system is defined by iterating the 
following steps: 

1. Select at random an agent i 6 a and a agent j G /?. 

2. Select the source of interaction: with probability B, 
agent i € a interacts with field Mp and agent j G 
f3 interacts with field M a , while with probability 
1 — B, i interacts with k G a selected at random 
and j interacts with I G f3 also selected at random. 

3. Calculate the overlap (number of shared compo- 
nents) between agent i G a and its source of interac- 
tion, given by d a = S/=i S xUi),y f,wheie yf = M l 
if the source is the field Mp, or y* = x^(k) if the 
source is agent k G a. Similarly, calculate the over- 
lap dp = E/=i^(j),<,/: where y f = M s a if the 

source is the field M a , or yf — Xp(l) if the source is 
agent / G /?. Here we employ the delta Kronccker 
function, 5 XtV = 1, if x = y; 5 x , y — 0, if x ^ y. 

4. If < d a < F, with probability %■ choose g such 
that x a a (i) ^ y 9 and set x 9 a (i) = y 9 ; if d a = 
or d a = F, the state x a (i) does not change. If 
< dp < F, with probability ^ choose h such 
that Xp(j) ^ y h and set Xp(j) = y h ; if dp = or 
dp = F, the state Xp(j) does not change. 

5. If the source of interaction is Mp[ a ], update the 
fields M a and Mp. 

The strength of each field M a and Mp is represented 
by the parameter B G [0, 1] that measures the probability 
for the agent-field interactions. 

Under the mutual coupling, both populations, a and 
(3 form domains of different sizes in the asymptotic state. 
A domain is a set of connected agents that share the 
same state. A homogeneous or ordered phase in a pop- 
ulation corresponds to d(i,j) — F, Vi,j. There are q F 
equivalent configurations for this ordered phase. In an 
inhomogeneous or disordered phase in a population sev- 
eral domains coexist. The sizes of these domains within 
each population are ranked by the index r: r = 1 corre- 
sponding to the largest domain, r = 2 indicates the sec- 
ond largest domain, etc. To characterize the collective 
behavior of the system, we define the following macro- 
scopic quantities: (i) the average normalized size (divided 
by N a \p]) of the domain in a[/3] whose size has rank r, 
denoted by (h) the probability that the largest 

domain in a[f3] has a state equal to Mp\ a ], designed by 

Figure [2] shows various of these quantities as func- 
tions of the parameter q, for different values of B. In 
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the absence of global fields (Fig. HJa)), i.e. B = 0, we 
have two uncoupled and independent subsets; each subset 
spontaneously reaches an ordered phase, characterized by 
S} t = 1 and Si = 1, for values q < q c , and a disordered 
phase, corresponding to ~ and Sp ~ 0, for q > q c , 
where q c is a critical point that depends on the subset 
size in each case, q c {a) ~ V Q [ill- Figure EJa) shows the 
asymptotic pattern in this case. 

For B — > and q < q c , each population reaches an or- 
dered state with S^ = 1, as shown in Fig.|^b) . However, 
in this situation the spontaneous order emerging in sub- 
set a for parameter values q < q c due to the agent-agent 
interactions competes with the order being imposed by 
the applied global field Mp. For some realizations of ini- 
tial conditions, the global field Mp imposes its state on 
subset a and, correspondingly, the field M a imposes its 
state on subset f3. As a consequence, both subsets reach 
the same state with M a = Mp. An asymptotic state cor- 
responding to this situation is displayed in Fig. [31(b) . We 
refer to this state as phase I. However, the ordered state 
in subset a[j3] does not always correspond to the state of 
the global field Mp\ a ] being applied to a[/3]. This is re- 
vealed by the probability P^(Mp) shown in Fig. [2Kb) that 
measures the fraction of realizations that the largest do- 
main in a has a state equal to Mp. We find P^(Mp) < 1 
for a range of values q < q c . Thus, in this case there is a 
probability that subsets a and /3 can reach ordered states 
different from each other, i. c, A lp ^ M a . Figure [3Jc) 
illustrates the asymptotic states in this case. We denote 
this situation as phase II. 

Figures [2Jc) and HJd) show both and Si as func- 
tions of q for greater values of B. The quantity in 
Fig.[2jc) displays a local minimum at some value of q that 
depends on B. This local minimum of S* is associated 
to a local maximum value of Si, such that + Si ss 1 
for q < q c . Therefore, two majority domains form in 
subset a for q < q c . Fig. [5Jc) also shows that the prob- 
ability P^(Mp) = 1, indicating that the state of the 
largest group in a is always equal to that imposed by 
the field Mp. But the second largest group that occu- 
pies almost the rest of subset a reaches a state different 
from Mp. Thus, the value of q < q c for which S& has a 
local minimum is related to the emergence of a second 
largest domain ordered against the global field Mp. The 
corresponding asymptotic pattern is shown in Fig. (3Jd). 
We call this configuration phase IV. Figure [U(d) reveals 
that, for larger values of B, various local minima of 
can occur at some values of q. This local minima of S} t 
correspond to local maxima of Si and to the emergence 
of a second largest domain in a ordered against the field 
Mp. The raise of a largest minority group at some values 
of q is a manifestation of the tendency towards the spon- 
taneous order related to the agent-agent interactions. For 
values q > q c . both populations reach disordered states 
VB, characterized by — Sp ~ 0. The disordered be- 
havior of the system is denoted by phase III and the 
corresponding pattern is displayed in Fig.[3](e). 

To characterize phase II, we plot in Fig.Uthc quantity 




FIG. 2: Si, Si, and P^{Mp) as functions of q, for differ- 
ent values of B. System size is N — 800 with partition 
N a = 0.67V. Each data point is the result of averaging over 
100 random realizations of initial conditions, (a) 5^ (open 
circles), Sa (solid circles); with B = 0. (b) Left vertical 
axis: S& (open circles); right vertical axis: Pa(Mp) (crosses); 
fixed B = 0.001. Phases I and II. (c) Left vertical axis: S& 
(open circles), (open squares); right vertical axis: P^{Mp) 
(crosses); fixed B — 0.05. Phases I and IV. (d) Left vertical 
axis: S^ (open circles), S„ (open squares); right vertical axis: 
P*(Mp) (crosses); fixed B = 0.25. Phase III occurs for values 
q > q c — 2500, independent of B. The left and the vertical 
axes in (b), (c), and (d) have the same scale. The bars in 
(c) and (d) indicate the probability <f> of finding a localized 
ordered state in the system as a function of q for the given 
value of the intensity B. 



er Q = (1 — P^(Mp))S^ that measures the average size 
of the largest domain in a whose state is different from 
Mp, as a function of q, for a fixed value B = 0.0005. 
For q <C q* ~ 10, the state of the largest domain in a 
corresponds to the state of the field Mp, i.e. P^(Mp) = 1 
and = 1, indicating the presence of phase I, and thus 
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FIG. 3: Each panel displays an asymptotic state (vertical 
axis) of the agents in the interacting populations a (upper 
part) and ft (lower part) vs. time (horizontal axis), cor- 
responding to a different phase in the system. Each value 
of the state variable of an agent is represented by a differ- 
ent color. Population sizes are N a — 0.6N, Np — 0.6N, 
with N = 800. (a) B = 0,q = 80 (no coupling), (b) 
B = 0.001,g = 80 (phase I), (c) B = 0.001, q = 100 (phase 
II). (d)B = 0.05, q = 80 (phase IV). (e) B = 0.25, q = 2500 
(phase III), (f) B — 0.05, q = 80 (localized ordered state). 



a a = 0. For q* < q < q c , the largest domain in a no 
longer possesses the state of the field Mp but another 
state non- interacting with this field, i.e. P^{Mp) < 1 
and = 1, and therefore a a > 0, characterizing phase 
II. For q > q c , S„ ->• and a a = 0. 




FIG. 4: The quantity a a = (1 - P^(M^))Si as a function of 
q for a fixed value B — 0.0005. The critical values q* and q c , 
as well as the regions where phases I, II, and III occur, are 
indicated. System size is N — 800 with partition N a — 0.6N. 
Each data point is averaged over 100 realizations of initial 
conditions. 

We note that phase II occurs for small values of B, 
where the time scale for the agent-agent interaction dy- 
namics is smaller than the corresponding time scale for 



the agent-field dynamics. This means that the state of 
the global field does not vary much in comparison to 
the changes taking place in the states of the agents and, 
therefore, the global field behaves approximately as a 
fixed external field with little influence on the system. 
As a consequence the system can spontaneously order in 
a state different from that of the global field if q < q c 
is sufficiently large, giving rise to phase II. For increas- 
ing values of B, the updating of the global fields and 
the agent-agent dynamics have comparable time scales 
and, therefore, the state of the fields corresponds to that 
of the largest domain in each subset, yielding regions of 
both phase I and phase IV. 

The collective behavior of either of the two subsets cou- 
pled through their reciprocal global fields can be charac- 
terized by four phases on the space of parameters (B, q), 
as shown in Fig. [5] for subset a: (I) a homogeneous, or- 
dered phase, for which ~ 1 and P^{Mp) = 1; (II) an 
ordered phase in a state orthogonal to the applied global 
field, such that ~ 1 and P«(M^) < 1; (III) a disor- 
dered phase for q > q c , for which ~ 0; and (IV) a 
partially ordered phase, where > and + w 1, 
P^(Mp) = 1, characterized by the emergence of a second 
largest domain ordered in a state different from field Mp. 
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FIG. 5: Phase diagram of population ot on the space of pa- 
rameters (B,q). System size is N = 800 with partition 
N a — 0.67V. Each data point is averaged over 100 realizations 
of initial conditions. The color code represents the value of 
the normalized largest domain size Si, from black (S a = 0) to 
white (S a = 1). The regions where the different phases occur 
are labeled and separated by slashed lines: phase I (both pop- 
ulations share same homogeneous state); phase IV (partially 
ordered, emergence of second group); phase III (disordered), 
and phase II (each population in a different homogeneous 
state). Localized ordered states can occur in the transitions 
from phase IV to phase I. 

The phase diagram of Fig. [5] reveals that the interac- 
tion through reciprocal, evolving global fields can lead 
to nontrivial effects in certain cases. For example, for a 
fixed value q = 20, the global field can impose its state 
to the system (phase I) only for a range of intermediate 
values of the intensity B. 

We have checked the behavior of the system for dif- 
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ferent population sizes N a and Np. Figure [6] shows the 
quantity as a function of q/N a with fixed coupling _B, 
for different values of N a . We see that the critical point 
for the transition to phase III scales as q c ~ N a , as ex- 
pected [4l|, and that the qualitative collective behavior 
represented in the phase diagram of Fig. [5] is indepen- 
dent of the sizes of the partitions into two populations. 
Since N a oc N, the collective behavior of the system is 
also independent of the size N, and q c <~ N, according to 

Fig.m 
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FIG. 6: Normalized size of largest domain S a as a function 
of q/N a with fixed intensity B — 0.1, for different population 
sizes: N a = 0.5iV (squares); N a = 0.56iV (circles); N a — 
0.8JV (diamonds). System size is N = 800. 



realizations of initial conditions, for different values of q, 
and with fixed B = 0.05 corresponding to Fig. 2(c). Fig- 
ureJTJa) exhibits the probabilities p(S^) and p(Sp) when 
either subset is in phase I with q = 10, characterized 
by the presence of one large domain whose size is of the 
order of the system size S^[Sl] ~ 1, in agreement with 
Fig. 3(b). FigureJTJb) shows p(S^) and p(Sl) associated 
to phase IV (q = 70), where the size of the largest domain 
in either subset never reaches the system size due to the 
appearance of a second group, as displayed in Figs. 2(c) 
and 3(d). Figure [7f c) shows the probabilities p(S^) and 
p(Sp) for q = 90. In this case either subset can reach an 
ordered configuration, S^fSo] ^ 1, or an inhomogeneous 
state (S^[Sp] < u). This corresponds to the appearance 
of localized ordered states in the system. For q = 115, 
we find again a probability distribution typical of phase 
I. 
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III. LOCALIZED ORDERED STATES. 

In addition to phases I and II that display homoge- 
neous states for both subsets a and j3, there are configu- 
rations where homogeneous states can take place in only 
one subset, while the other is inhomogeneous, for some 
values of parameters. We refer to this configuration as 
localized ordered states. These states are characterized 
by Si [S}} = 1 and Sj [S^ = u < 1. Figure Eff) dis- 
plays the asymptotic state of the system in this case. In 
contrast to the four phases that can be characterized in 
a subset, the ordered collective states can only be de- 
fined by considering both subsets simultaneously, i.e., it 
requires the observation of the entire system. 

To elucidate the nature of these states, we calculate 
the probability <f> of finding a localized ordered state in 
the system as a function of q in Figs. 2(c) and 2(d), em- 
ploying the criterion u < 0.6. In both figures, there are 
ranges of the parameter q where localized ordered states 
can occur; the probability </> is maximum near the values 
of q that correspond to local minima of (and local 
maximum values of S^)- Figure [7] shows the probability 
distributions p(S^) and p(S^), Vr, of the normalized do- 
main sizes for both subsets a and /?, calculated over 100 



FIG. 7: Probability distributions p(S„) and p(S%), Vr, of nor- 
malized domain sizes for both populations ot (black bars) and 
P (grey bars), calculated over 100 realizations of initial con- 
ditions, with fixed intensity B = 0.05 and for different values 
of the number of options q. (a) q = 10 (phase I); (b) q = 70 
(phase IV); (c) q = 90 (localized ordered states); (d) q — 115 
(phase I). 



The localized ordered states are analogous to chimera 
states observed in two populations of dynamical oscilla- 
tors having global or long range interactions, where one 
population in a coherent state coexist with the other in a 
incoherent state [HI, 39] . In a chimera state, one part of 
a spatially extended system presents a coherent or syn- 
chronized behavior while another part is desynchronized 

MM. 



Note that the regions of parameters where localized or- 
dered states can emerge in our system lie between phase 
IV and phase I states. In fact, the configuration of local- 
ized ordered states shares features of both phase I and 
phase IV; they can be considered as transition configu- 
rations between phase IV and phase I states. 
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IV. DISCUSSION. 

We have investigated the collective behavior of a sys- 
tem consisting of two populations of social agents, mutu- 
ally coupled through global fields, as a model for cross- 
cultural interactions via mass media. Specifically, we 
have employed Axelrod's model for social influence as 
the interaction dynamics. 

The global interaction field associated to each popu- 
lation corresponds to the statistical mode of the states 
of the agents. In the context of social dynamics, this 
global autonomous held can be interpreted as mass me- 
dia messages about "trends" or stereotypes originated in 
one population that are transmitted to the other pop- 
ulation. Thus, our system can represent cross cultural 
interactions between two separated social groups, each 
with its own internal dynamics, but getting information 
about each other solely through their mass media mes- 
sages [||. 

We have found several phases on either subset depend- 
ing on parameter values: two homogeneous phases, one 
having the state of the global field acting on that subset 
(phase I), and the other consisting of a state different 
from that reached by the applied global field (phase II); 
a partially ordered phase characterized by the emergence 
of a second largest domain ordered in a state different 
from the global field (phase IV); and a disordered phase 
(III). 

States similar to phases I, II, and III are also ob- 
served for some regions of parameters in a system of so- 
cial agents subject to an external fixed field [42|]. In the 
present model, for small values of B, the global evolving 
field varies very slowly in comparison to the changes in 
the states of the agents in a subset due to their mutual in- 
teractions. In this case, the global evolving field behaves 
as a fixed external field acting on the population. 

However, for larger values of B, the adaptive nature 
of the global fields induce two new phenomena in some 
range of q < q c on each population. One is the emergence 
of a largest minority group of appreciable size having a 
state different from that of the applied field (phase IV). 
The other corresponds to the appearance of localized or- 
dered states when the entire system is observed, consist- 
ing of one population in a homogeneous state and the 
other in an disordered state. These configurations occur 
with a probability that depend on both B and q and ap- 
pear as transitions states from phase IV to phase I. These 



localized ordered states arc analogous to the chimera 
states that have been found in networks of coupled oscil- 
lators having global interactions, where a subset of the 
system reaches a coherent state while another subset re- 
mains incoherent [H, |39| . The discovery of such chimera 
states has fundamental implications as it shows that lo- 
calized order and structured patterns can emerge from 
otherwise structureless system. As noted in Ref. [38j |. 
analogous symmetry breaking is observed in dolphins and 
other animals that have evolved to sleep with only half 
of their brain at a time: neurons exhibit synchronized 
activity in the sleeping hemisphere and desynchronized 
activity in the hemisphere that is awake [43| . 

From a social perspective, our model shows that cross 
cultural reciprocal interactions through mass media do 
not always lead to the imposition over one population 
of the cultural trends being transmitted by the media of 
another population. A group possessing a cultural state 
different from that of the mass media message can spon- 
taneously emerge in the first population. Under some cir- 
cumstances, such group can encompass the entire popu- 
lation (phase II), or it can constitute the largest minority 
in that population (phase IV). 

The behaviors reported here should also be expected in 
other non-equilibrium systems possessing non-interacting 
states, such as social and biological systems whose dy- 
namics usually possess a bound condition for interaction 
[3l| . This includes models of motile elements in popula- 
tion dynamics, such as swarms, fish schools, bird flocks 
and bacteria colonies [30l . |44T - |48| . Future extensions of 
this work involves the consideration of complex network 
structures within each population and the investigation 
of communities, where the interaction between popula- 
tions occurs through a few elements rather than a global 
field. 
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